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Abstract 

We study the production process e + e~ — > X1X2 anc ^ the subsequent 
decay X2 ~ * Xi^ + ^ with polarized e + and e _ beams, including the 
spin correlations between production and decay. We work out the 
advantages of polarizing both beams. We study in detail the angular 
distribution and the forward-backward asymmetry of the decay lepton 
as well as the opening angle distribution between the decay leptons. 
We investigate the dependence on the masses of e_L and £r and on 
the mixing character of the neutralinos. In particular we study the 
dependence on the gaugino mass parameter M\. 
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1 Introduction 



The search for supersymmetric (SUSY) particles and the determination of 
their properties will be one of the main goals of a future e + e~ linear collider. 
Particularly interesting will be the experimental study of the neutralinos, 
which are the quantum mechanical mixtures of the neutral gauginos and 
higgsinos, the SUSY partners of the neutral gauge and Higgs bosons. In the 
Minimal Supersymmetric Standard Model (MSSM) there are four neutrali- 
nos Xii i = 1, ••• ,4. In extensions of the MSSM there may be more than 
four neutralinos. Usually the lightest neutralino Xi i s the lightest SUSY 
particle LSP. Therefore, the production of the lightest and the second light- 
est neutralino can presumably be studied at an e + e~ linear collider with 
a CMS energy = 500 GeV. The aim will be to precisely determine the 
SUSY parameters of the neutralino system. By a detailed study of the neu- 
tralino system one can also examine the question whether the MSSM or 
another SUSY model is realized in nature. Models with an extended neu- 
tralino sector have been discussed in Q . Recently a method for determining 
the SUSY parameter M 2 , \x and tan f3 by measuring suitable observables in 
chargino production e + e~ — > xfxj (hj = U2) has been proposed in ||. 
The process e + e~ — ► xtxl " ■> X.t ~^ Xk^ +U including the full spin correlation 
has been studied in ||]. 

In previous papers (see || and references therein) the cross sections 
of e + e~ — ► XiXp hj = 1 • • • 4, and the branching ratios and energy dis- 
tributions of neutralino decays were studied, which do not depend on spin 
correlations. In the calculation of the decay angular distributions, however, 
one has to take into account the spin correlations between production and 
decay of the neutralinos. In Q we have studied the process e + e~ — > XiXp 
i,j = 1...4, with unpolarized beams and the subsequent leptonic decays 
Xi ~^ Xk^ + f-~ including the complete spin correlations. In j/J we have given 
the complete analytical formulae for longitudinally polarized beams, fully 
including the spin correlations between production and decay. The formu- 
lae have been given in the laboratory system in terms of the basic kinematic 
variables. 

In the present paper we extend our analysis in ]7| and study neutralino 
production in the case that both beams are polarized. We show that by a 
suitable choice of the e + beam polarization not only higher cross sections but 
also additional information can be gained. We give numerical predictions 
for the process e + e~ — > X1X2 with X2 ~^ Xi e+e ~ including the full spin 
correlations between production and decay. The framework of our studies 
is the MSSM. We first assume the GUT relation M 1 /M 2 = § tan 2 @ w for 
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the gaugino mass parameters (note that in Refs. [Q, ^, [?J we used the 
notation M' and M for M\ and M2). We consider two gaugino-like and one 
higgsino-like scenario and study the dependence of the cross section, the 
forward-backward asymmetry of the decay electron, and the opening angle 
distribution of the decay lepton pair e + and e~ on the beam polarizations 
and on the masses of the exchanged and en- Then we relax the GUT 
relation between M\ and M2 and study the M\ dependence of a(e + e~ — ► 
X1X2) x BR(x2 ~^ Xi e+e ~) an d the forward-backward asymmetry of the 
decay electron for various beam polarizations and slepton masses. We also 
discuss which observables are suitable for the determination of the neutralino 
parameters and selectron masses. 



2 Spin correlations between production and decay 

Both the helicity amplitudes Xp Xj for the production process 
e~(pi)e + (j>2) -» X°(^3, Aj)x°(p4, Aj) and the helicity amplitudes T D;At 
and Tjj\. for the decay processes Xi(.P3,^i) ~^ Xk(PsK + \p&)^~ (Pr) an d 
X°(p4,Aj) — > X?(P8)^ + (P9)^~(Pio)> respectively, receive contributions from 
Z° exchange in the direct channel and from ez, r exchange in the crossed 
channels. 

The amplitude squared of the combined process of production and decay 

is: 

\ T \ 2 = |A(x°)| 2 |A(x°)| 2 pp A ' iXj p DtX '. x .p DiX ', x . (summed over helicities).(l) 
It is composed of the (unnormalized) spin density production matrix 

A/AiA A A A A A * , , 

Pp l 1 1 3 =Tp lX] Tp 3 , (2) 
the decay matrices 

P D ,X'.\ = Td ^ T *D A' alld P D,X'.Xj = T D,h T * D A' ' ^ 
and the propagators 

A(Xi,i) = Vbi,4 - ml 3 + im id T ijjt ]. (4) 

Here P34, Ajj, my and Tij denote the four-momentum squared, helicity, 
mass and total width of Xij- For these propagators we use the narrow-width 
approximation . 
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Introducing a suitable set of polarization vectors for each of the neutrali- 
nos the density matrices can be expanded in terms of Pauli matrices a a : 

+ V E <^pm + E ^a<.a' s ^°x°))> ( 5 ) 

PD,\'\i = ^% D (x°)+Y,°i%vum (6) 

a 

^,a'a, = (Va^^^ + E^a'a,^^ )) ^,b=l,2,3. (7) 

3 J J b ^ 

We choose the polarization vectors such that Sp(x^) describes the trans- 
verse polarization in the production plane, Sp(x° •) denotes the polarization 
perpendicular to the production plane and Sp(x^) describes the longitudi- 
nal polarization of the respective neutralino. Sp (xi xj) is due to correlations 
between the polarizations of both neutralinos. The complete analytical ex- 
pressions for the production density matrix and for the decay matrices are 
given in (?]]. 

If CP is conserved the neutralino couplings are real. It can be seen in ]7| 
that in this case, due to the Majorana character of the neutralinos, the quan- 
tities P, Sp, Sp 1 , Sp , Sp , Sp are forward-backward symmetric, whereas 
the quantities Sp, Sp, Sp , Sp are forward-backward antisymmetric. 

The amplitude squared of the combined process of production and decay 
can be written as: 

|T| 2 = 4|A(x°)| 2 |A(x°)| 2 (P(x°X?) J D(x°) J D(x J °) + E ^ a p(x^ a D (xf)D(x° 3 ) 

a=l 

+ ESp(x°)S^(x°) J D(x°) + E ^f(x^Ux^ b D (x°,))- (8) 

6=1 a,b=l 

The differential cross section in the laboratory system is then given by 
da e = J-|T| 2 (2^) 4 5 4 ( Pl +p 2 -J2Pi)dtips(j>3 ■ ■ -Pio), (9) 

i 

<ilips(p3, . . . ,Pio) is the Lorentz invariant phase space element. 

If one neglects all spin correlations between production and decay only 
the first term in (||) contributes. The second and third term in (||) describe 
the spin correlations between the production and the decay process and the 
last term is due to spin-spin correlations between both decaying neutralinos. 
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3 Numerical analysis and discussion 



In the MSSM || the masses and couplings of neutralinos are determined by 
the parameters M\, M 2 , fi, tan/3, which can be chosen real if CP violation 
is neglected. Moreover, one usually makes use of the GUT relation 

Mi = ^M 2 tan 2 Q w . (10) 

The neutralino mass mixing matrix in the convention used can be found in 
§- 

The total cross section of e + e — ► XiXj an d the decay rate Xi ~ * Xk e+e 
further depend on the masses of and §r. In the following numerical 
analysis we study neutralino production and decay in three scenarios, which 
we denote by Al, A2, and B. The corresponding parameters are given in 
Table [Tj. 

Scenario Al corresponds to that in ||. In this scenario Xi is P - like and 
X2 is PF 3 -like. The selectron masses are = 176 GeV, = 132 GeV. 
Scenario A2 differs from Al only by the mass of ex- In scenario B the 
same masses as in scenario Al are taken for x?> X2> Xi an d X2 are ; 

however, higgsino-like due to the choice fi < M 2 . 

3.1 Effects of beam polarizations on the total cross section 

The cross section cr(e + e _ — ► X1X2) is shown in Figs, la, lb, and lc as a 
function of the longitudinal beam polarizations P 3 for electrons and P 3 . for 
positrons, for scenario Al, A2 and B, respectively, (with P± = {—1,0, 1} for 
{left-, un-, right-} polarized ). The cross section a(e + e~ —> X1X2) is shown 
at \fs = (m^o + m^o) + 30 GeV. 

v v Al X-2' 

In Figs, la, lb, and lc the white area is covered by an electron polar- 
ization |P 3 | < 85% and a positron polarization \Pl\ < 60%. The cross 
section can be enhanced by a factor 2-3 by polarizing both beams. Theo- 
retically, for pure gaugino-like neutralinos and mg L 3> mg R {rrie L <C m e -„) 
and P 3 = +1, P| = —1 (P 3 = —1, P 3 . = +1), the cross section could be 

enlarged by a factor 4. For pure higgsino-like neutralinos and P£ = +1, 
p3 = -1 (p3 _ _i P 3 = +1 ) this factor would be 17 (2.3) gg. 

One clearly recognizes, Figs, la, lb, and lc, the sensitive dependence 
of the cross section on the selectron masses mg L and as well as on 
the mixing character of Xi an d X 2 - in scenario A2 with mg L 3> rrie R an d 
gaugino-like x\ 2 i one expects the largest cross section for P 3 = +1 and 
P? = —1, see Fig. lb. In scenario B the cross section is governed by Z° 
exchange and is therefore rather symmetric for P 3 = ± «-» P 3 . = =p. The 
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beam polarizations are a useful tool for getting more information about m. 



and ingu- 



lf the polarizations of both beams are varied, the relative size of the 
cross sections strongly depends on the mixing character of both neutralinos 
X?2 an d on the selectron masses [[R]]. If x± an d \\ are pure higgsinos, one 
obtains for |P£| = 85% and |P?| = 60% the sequence 



o- + > at' > of > a™ > a+° > *T > (11) 
Here a e = a(e~e + — >■ X1X2) x BR(x2 — ► an d ( — I - ) etc denotes 



the sign of the electron polarization P£ and of the positron polarization P? 



respectively. 

If Xi an d X2 are pure gauginos the order of the cross sections depends 
on the relative magnitude of the selectron masses and mg R ■ For mg L ^> 
rrie R only right selectron exchange contributes, and one obtains 

of > > °e° > °t + > °~e~ > V7° > °e + , (12) 

whereas for 3> mg L , one gets: 

> *e° > ^° > <?r > °i + > °-t° > °t~- (is) 

The case of a heavy right slepton may be realized in extended SUSY models 



([11] and references therein). 

Comparing dill) and ([l3]) shows that polarizing both beams allows one 
to distinguish between a higgsino-like scenario and a gaugino-like scenario 
with dominating exchange. This is not possible if only the electron beam 
is polarized. 

In Table ^ we show the cross sections for various polarization configu- 
rations for our scenarios Al, A2, and B at = (rn^o + m^o) + 30 GeV 
for P£ = 0, ±85% and P+ = 0, ±60%. For the higgsino-like scenario B the 
sequence of the cross sections coincides with ( |TT1) . One notices that one ob- 
tains the same ordering of the polarized cross sections for the gaugino-like 
scenario Al. This shows that the relative size of the cross sections sensi- 
tively depends on the mass difference between ei and Sr, which is rather 
small in scenario Al, — mg H = 44 GeV. Comparing the sequence of 
cross sections for the gaugino-like scenario A2, see Table |2[ with that for 
pure e.R exchange, see (|i~2j), one sees a small influence of ei exchange despite 
the rather high ez, mass, = 500 GeV. 
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M 2 


V- 


m- eL 


m e R 


Xi 


X9. 


T x°, 


u 12 


feifk 


III J 12 


Al 


152 


316 


176 


132 


71 


130 


25E-6 


-.02 


-.20 


-.12 


A2 


152 


316 


500 


132 


71 


130 


15E-6 


-.02 


-.20 


-.12 


B 


250 


125 


176 


132 


71 


130 


369E-6 


+.39 


2E-5 


.027 



Table 1: Parameters, masses, and total X2 width (in GeV) and couplings in 
scenarios Al, A2, and B for tan/3 = 3. 



3.2 Lepton forward— backward asymmetry 

Owing to the Major ana character of the neutralinos the angular distribu- 



tion of the production process is forward-backward symmetric [12]. The 
angular distribution of the decay lepton, however, depends sensitively on 
the polarization of Xi- Since the longitudinal polarization £p and the trans- 
verse polarization £p of jn are forward-backward antisymmetric, the lepton 
forward-backward asymmetry Afb of the decay lepton may become quite 
large. The lepton forward-backward asymmetry Afb is defined as 

_ cx e (cose e > 0) - o- e (cos9 e < 0) 
FB ~ a e (cos G e > 0) + o- e (cos 9 e < 0) ' ( ' 

where a e is a short-hand notation for a e = a(e + e~ — > X1X2) x BR(x2 ~^ 
Xie + e~). We will show Afb not too far from threshold because it decreases 
with y/s for fixed neutralino masses. 

In Figs. 2a, 2b, and 2c we show Afb of the decay electron as a function 
of the electron and positron polarizations for the scenarios Al, A2, and B, 
respectively, at yfs = (m^o+m^o)+30 GeV. First one notices that polarizing 
suitably both beams gives a larger asymmetry. Actually, in the scenarios 
Al and B Afb turns out to be practically zero for both beams unpolarized. 
The figures also exhibit a very different pattern. When comparing Fig. 2a 
with Fig. 2b, the different behaviour is due to the different masses of e^. 
In the higgsino scenario B the asymmetries are much smaller. Notice again 
the symmetry of P£ = ± <-> Pi = =p in scenario B as already observed in 
the total cross section. Measuring the lepton forward-backward asymmetry 
Afb i n addition to the total cross section strongly constrains the selectron 
masses mg L and mg R as well as the mixing properties of x± and x% 

We also studied the dependence of the lepton forward-backward asym- 
metry Afb on 1/s. For y/s 3> (m^o + m^o) the angular distribution of the 



decay lepton is essentially the same as that of the decaying neutralino X2 
Therefore the lepton forward-backward asymmetry practically vanishes. 
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ffe/fb 


9.5 


6.0 
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(-+) 


(+-) 


(-0) 


(00) 


(+0) 


(— ) 


(++) 


ff e /fb 


21.7 


19.0 


14.3 


13.5 


12.7 


6.8 


6.4 



Table 2: Polarized cross sections a e = a{e + e — > X1X2) x BR(x2 ~^ 
Xie + e~)/fb at ^fs = m^o+m^o + 30 GeV in scenarios Al, A2, and B, see Ta- 
ble |l|, for unpolarized beams (00), only electron beam polarized (—0), (+0) 
with Pf = ±85% and both beams polarized with Pf = -85%, P 2 3 = +60% 
(-+) and Pf = +85%, P| = -60% (+-). 



3.3 Opening angle distribution 

The opening angle distribution between the two leptons from the decay of 
one of the neutralinos, X2 ~^ ^ + ^~Xi> is independent of the spin correlations 



due to the Majorana nature of the neutralinos [14]. Therefore, it factorizes 
into the contributions from production and decay. For the same reason this 
is also valid for the energy distribution of the neutralino decay products. 
For both distributions it is suitable to parametrize the phase space by the 
scattering angle between the incoming e~(p\) beam and the outgoing neu- 
tralino x^iPi)^ the azimuthal angle ^^o e - between the scattering plane and 

the (x2^ _ ) _ plane and the opening angle 0_| between the leptons £ + and l~ 

from the decay of the neutralino x% Since the phase space is independent of 
the azimuthal angle, the contributions of the transverse polarizations of the 
neutralino vanish after integration over &%og- . As for the longitudinal polar- 
ization Tp, the Majorana character of the neutralino is crucial. If CP is con- 
served Tp is forward-backward antisymmetric, Tp(— cos0) = — Ep(cosO), 
so that the contribution of the longitudinal polarization vanishes after inte- 
gration over the scattering angle [ fH|] . 

In Figs. ||a and ||b we show the distribution of the angle 0^ between the 

decay leptons from e + e~ — > X1X2) X2 ~^ e+e ~Xi i n t ne scenarios Al, A2 and 
B at = (m^o + Tn^o) + 30 GeV and for various beam polarizations. One 

notices that the shape of the 0_| distribution is mainly determined by the 

mixing character of the neutralinos. The selectron masses mainly influence 
the size of the cross section. Due to the factorization of production and decay 
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the beam polarizations have no influence on the shape of this distribution. 



3.4 Dependence on Mi 

So far we have used the GUT relation (10) for the gaugino masses. In the 



following we will be more general and not use this relation |||, [l3|, 15, 
We will discuss the Ml dependence of the cross section and the forward- 
backward asymmetry of the decay electron All other parameters are 
chosen as in scenario Al except the mass of e/j, = 161 GeV. The neu- 
tralino masses as well as the Z^x^X^j couplings and the jQil couplings 
f^ ,R depend on M\. 

In Fig. | we show the neutralino masses as function of M\. The grey 
areas are excluded by the constraints m^o < m~±, m^o > 35 GeV. We 

J Xi Xi ' Xi 

see that in the interval between —130 GeV< Mi < M2 m^o depends very 
strongly on Mi, whereas all other neutralino masses are nearly independent 
of Mi. On the other hand, in the region M2 < Mi < \u\ only m^o depends 

A 2 

strongly on Mi. 

In the formulae for the cross section a{e + e~ —* X1X2) aim f° r the decay 
X2 ~^ Xi^ + £~ the products of couplings jf-yjf 2 and f R f$> enter. We therefore 
show for these products the dependence on Mi in Fig. [5|. We do not consider 
values I Mi I > 160 GeV, where m^o > mg R , because then the two-body 
decay \\ ~~ * + e would be possible. One observes a strong variation 
of f R f^ and f({fi2 for M\ > 80 GeV. In particular, f R f^ has a positive 
maximum at 140 GeV, whereas ff\ff 2 ls zero a t M\ = 120 GeV and reaches 
large negative values for Mi > 160 GeV. We therefore have the following 
regions: f e L j[ 2 > f R f R > for -200 GeV< Mi < 80 GeV, \f R f R \ > tfjg 
for 110 GeV< Mi < 140 GeV, and l/^/ijI > f&fi for M i > 150 GeV - The 
coupling O12 is small in this gaugino-like scenario. 



Fig. |6|a exhibits the Mi dependence of cr(e + e — > X1X2) x BR(x2 
o0\ 



e ■ e~Xi) at ^/s = (m^o + rn^o) + 30 GeV in the region 40 GeV< 
Mi <160 GeV for various beam polarizations. Since the masses of ez, and 
eR are in this case comparable, the curves reflect the behaviour of fkfh an d 
f R f(2 °f Fig. |5[ A left (right) beam polarization of the electron (positron) 
selects the exchange, while the maximum in the curve with a right (left) 
electron (positron) polarization is due the maximum of f R f(2 m Fig- H- 

Fig. [6|b shows the analogous curves for a heavy ei (,irie L = 500 GeV) and 
all other parameters as in Fig. |||a. One clearly sees that the &l exchange is 
strongly suppressed, and one obtains higher cross sections for right polarized 
e~ beams. 
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We have also studied the M\ dependence of the forward-backward asym- 
metry Aps of the decay electron, eq. (14). It is shown in Fig. f?]a for 
rrie L = 176 GeV, and in Fig. for mg^ = 500 GeV. One notices a strong 
variation with M\ and a strong dependence on the beam polarizations. 
Comparing Fig. |?]a and Fig. |7|b, one observes a very pronounced difference 
of the forward-backward asymmetry of the decay electron in the region 
40 GeV< Mi <100 GeV. This is due to the suppression of exchange in 
Fig. |7|b. The beam polarizations enhance the effect considerably. The peak 
at Mi 120 GeV is again due to the maximum of f^f^- 



4 Conclusions 

The objective of this paper has been twofold. Firstly, we have studied the 
advantage of having both the e~ and the e + beam polarized. If the polariza- 
tions of e~ and e + are varied, the relative size of the cross sections depends 
significantly on the mixing character of the neutralinos and on the masses of 
&l and eR. By an appropriate choice of the polarizations one can obtain up 
to three times larger cross sections than in the unpolarized case. Secondly, by 
taking into account the full spin correlations between production and decay, 
we have studied the angular distribution, as well as the forward-backward 
asymmetry of the decay electron e + e~ — ► XiX2> X2 ~~ > Xi e+e ~- Measur- 
ing this asymmetry for various beam polarizations strongly constrains the 
masses of and en and the mixing properties of the neutralinos. We have 
also studied the dependence on the gaugino mass parameter M\. For a 
determination of Mi the use of polarized e + and e~ beams would be very 
helpful. Due to the Majorana character of the neutralinos the opening angle 
distribution between the decay leptons is independent of spin correlations. 
It is very sensitive to the mixing character of the neutralinos, whereas its 
shape is only weakly dependent on the selectron masses. 
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Figure 1: Contour lines of cross sections 
cr(e + e" -> XiX^/fb at = (™x? + 
m^o) + 30 GeV in a) scenario Al with 
mel = 176 GeV, mg R = 132 GeV, in 
b) scenario A2 with m^ L — 500 GeV, 
mg R = 132 GeV, and in c) scenario B with 
m gi = 176 GeV, m~ eii = 132 GeV. The lon- 
gitudinal beam polarization for electrons 
(positrons) is denoted by P? (Pi). The 
shaded region is for IP 3 ] > 85%, |P+| > 
60% (dashed-linc if only electron beam po- 
larized) . 
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Figure 2: Contour lines of the forward- 
backward asymmetry of the decay electron 
A FB /% of e+e- - xlxlxl - xle + e~ 
at y/s = (m^o + m%a) + 30 GeV in a) 
scenario Al (mg £ = 176 GeV, mg R = 
132 GeV), in b) scenario A2 (mj t = 
500 GeV, mt R = 132 GeV), and in c) 
scenario B (mg^ = 176 GeV, rri£ L = 
176 GeV). The longitudinal beam polar- 
ization for electrons (positrons) is denoted 
by (P+)- The shaded region is for 
\Pl\ > 85%, \Pl\ > 60% (dashed-line if 
only electron beam polarized). 
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Figure 3: Opening angle distribution of the decay leptons from e + e~ — ► 
XiX2; X2 ~ * e+e ~Xi m a ) scenario Al (mg^ = 132 GeV, mg L = 176 GeV), A2 
(ra^ = 132 GeV, mg L = 500 GeV) and in b) scenario B (ra^ = 132 GeV, 
rrie L = 176 GeV) at = (m^o + rn^a) + 30 GeV for unpolarized beams 



(00), for P_ = -85%, P 4 
(H — ), respectively. 



+60% (-+) and for P_ = +85%, P 4 
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GUT: 78.7 Mi /GeV 

Figure 4: Neutralino mass spectrum as a function of M\ in the range 
—500 GeV< Mi < 500 GeV. The grey areas are excluded by the constraints 
m^o < m^t, Tn^o > 35 GeV. All other MSSM parameters as in scenario Al; 
GUT relation © gives M x = 78.7 GeV. 
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Figure 5: Slepton couplings fj{fg2 and fHf^ as a function of Mi in the range 
-200 GeV< Mi < 200 GeV. All other MSSM parameters as in scenario Al; 
GUT relation @ gives M 1 = 78.7 GeV. 

cr e /fb a) cr e /fb b) 




40 60 80 100 120 140 160 40 60 80 100 120 140 160 

Mi / GeV Afi /GeV 

Figure 6: Cross sections a e = a(e + e~ — > X1X2) X BR(X2 ~~ * Xi e+e ~) at yfs = 
m^o + m^o + 30 GeV as function of gaugino parameter Mi for unpolarized 

beams (00), for only electron beam polarized (—0), (+0) with = ±85% 
and for both beams polarized ( — h), (H — ) with P_ = =f85%, P + = ±60%. 
The slepton masses are a) mg L = 176 GeV, nie R = 161 GeV, and b) mg L = 
500 GeV, rrig R = 161 GeV; the other SUSY parameters as in scenario Al. 
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Figure 7: Forward-backward asymmetry of decay electron Afb of e + e~ — > 
XiX2> X2 — * Xi e+e ~ at \/s = n-^o + m^o + 30 GeV as function of gaugino 
parameter Mi for unpolarized beams (00), for only electron beam polarized 
(-0), (+0) with P 3 = ±85% and for both beams polarized (— +), (H — ) 



with P_ = =f85%, P + = ±60%. Slepton masses are a) 



in, 



l-'L 



176 GeV, 



nie R = 161 GeV, and b) mg L = 
parameters as in scenario Al. 



500 GeV, mg R = 161 GeV; the other SUSY 
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